Khinchin theorem and anomalous diffusion 
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A recent paper [M. H. Lee, Phys. Rev. Lett. 98, 190601 (2007)] has called attention to the fact 
that irreversibility is a broader concept than ergodicity, and that therefore the Khinchin theorem 
[A. I. Khinchin, Mathematical Foundations of Statistical Mechanics (Dover, New York) 1949] may 
fail in some systems. In this Letter we show that for all ranges of normal and anomalous diffusion 
described by a Generalized Langevin Equation the Khinchin theorem holds. 



1. Introduction. — The ergodic hypothesis (EH) states 
that the time and ensemble averages of phase variables 
exist and are equal for a stationary system. This hy- 
pothesis enables us to calculate thermodynamic quanti- 
ties from the equations of motion of the particles and 
is crucial for the proof of basic theorems in statistical 
mechanics [H, Q, OJ 0, HI ■ One of those theorems is the 
Khinchin theorem (KT) 2], of great importance, since it 
relates the ergodicity of a variable p to the irreversibility 
of its autocorrelation function. In a recent work [l[ based 
on the well established method of recurrence relations Q , 
it was asserted that, contrary to the KT, irreversibility is 
not a sufficient condition for ergodicity; in other words, 
the KT may not be valid for all systems. This demon- 
stration poses a new challenge, i.e., that of determining 
in which systems the KT is valid. 

Most of the experimental situations in which the EH 
does not hold arise in complex nonlinear or far from equi- 
librium structures where detailed balance is not fulfilled. 
A few examples are found in supercooled liquids 0, Q , 
glasses @, m, E3] and blinking nanocrystals [ll|. The 
majority of those systems, however, apparently do not 
have an easy analytical solution, with an exception hav- 



ing recently been reported in 12j. On the other hand, 
even anomalous diffusion treated with the Generalized 
Langevin Equation (GLE) can present closed solutions 
for the main expectation values, and can be used as a sim- 
ple laboratory for the discussion of those properties. As 
we shall see, in this case, we can give a full description for 
the validity of the KT, even when the EH breaks down. 
Although anomalous diffusive processes are present also 
in the context of deterministic Hamiltonian maps [l3j |. 
we shall focus our attention on the GLE formalism. 

In this Letter we work with stochastic processes and 
show that the KT is valid for all ranges of anomalous 
diffusion described by a GLE, even if this condition fails 
for Hermitian systems [ij. 

2. Khinchin's theorem. — Let p be a dynamical stochas- 
tic variable (e.g., momentum operator) of a classical par- 



ticle. The relaxation function can be written as 

o P (t) 



R(t) = 



C p (0) 



(1) 



where the autocorrelation function of p is given by 
C p {t) = (p(t)p(O)) - (p(t))(p(0)), and (...) stands for 
an ensemble average. Explicitly, the KT states that if 



R(t -> oo) = 0, 



(2) 



then p is ergodic 0]. In other words, the irreversibility 
condition is a necessary and sufficient condition for the 
validity of the EH. In Ref. it was claimed that for a 
system to be ergodic, it is necessary that < W < oo, 
where 



W 



R(t')dt'. 



(3) 



If this condition is true and if p refers to a Hermitian 
system, then irreversibility is not a sufficient condition 
for ergodicity, since the validity of Eq. |J5J| does not imply 
that W < oo. 

3. Diffusion phenomena. — Diffusive dynamics is usu- 
ally analyzed using the mean square displacement of the 
particles, which behaves in general as 

(lx(t) - (x(t))} 2 ) oc r, 

where the exponent a classifies the different types of dif- 
fusion: subdiffusion for < a < 1, normal diffusion for 
a = 1, and superdiffusion for 1 < a < 2; for a = 2 the 
process is called ballistic IBI. U3.[l5l lij 17 1. According to 
Kubo's linear response theory [181 ]. the diffusion constant 
is given by 



D= lim ±(l x (t)-(x(t))} 2 ) = ^r 



R(t')dt', 



(4) 

where m is the mass of the particle. Thus, for normal 
diffusion < D < oo, for subdiffusion D — 0, and for 
superdiffusion D = oo. According to the result of [l[, 
ergodicity is only ensured for normal diffusion. Since at 



2 



the time the KT was formulated most of the processes 
studied displayed normal diffusion with exponential re- 
laxation, it is quite natural to inquire if the KT will be 
affected, for example, in the slow relaxation dynamics 
that occurs in anomalous diffusion. In the anomalous 
regime, power laws, stretched exponentials, and Bessel 
functions are only a few examples of the vast functional 
behavior that is possible for relaxation [HI, . 

A general description of the diffusion dynamics can 
be given by means of the GLE, which was developed by 
Mori [2(j | using a Hermitian formulation, allowing to de- 
scribe all diffusive regimes including those beyond the 
Brownian limit. The GLE for a single particle in the 
absence of a net external force can be written as 



dpjt) 
dt 



= - [ n(t - t') P (t')dt' 

Jo 



(5) 



where LT(t) is the memory function, and r/(t) is a random 
force of zero mean. Besides this, the noise is uncorrelated 
with the initial value p(0), (rj(t)p(0)) = 0, and obeys the 
fluctuation-dissipation theorem (FDT) [181 ]: 



W*MO> 



where (. . .) eq is an average over an ensemble in thermal 
equilibrium. The solution of the GLE is 



p(t) = p(0)R(t) + R{t - t')r](t')dt' , (6) 
Jo 

where R(t) can be obtained by multiplying Eq. ([5|) by 
p(0) and taking the ensemble average, 



dR{t) 
dt 



R(t')U(t - t')dt' 



(7) 



whose Laplace transform yields R(z) = l/[z + n(z)]. In 
this sense, x(*) = —dR(t)/dt acts as a response func- 



tion 



18] 



4- Ensemble average and equilibrium condition. — If a 
system is ergodic and there are no external forces, ther- 
mal equilibrium should be observed in a time t 3> T c , 
where r c is a relaxation time. Then, the distribution 
function of p approaches the equilibrium distribution in 
the limit t — > 00, and the mean energy converges to the 
equilibrium value, (p 2 (t — > 00)) = (p 2 ) eq . 

For any initial distribution of values, p(0), it is possible 
to obtain the temporal evolution of the moments (p n (t)), 
with n = 1, 2, . . .. The first moment is obtained directly 
by taking the ensemble average of Eq. |6j): 



(p{t)) = (p(0))fl(i). 



(8) 



Taking the square of Eq. ([6]) and performing an ensemble 
average, we get 

(p 2 (t)) = (p 2 ) eq + R 2 (t) [{p 2 (0)) - . (9) 



Consequently, we see that the knowledge of R(t) allows 
one to describe completely these averages. Equations 
© and (JHJ) are sufficient to show the condition of equi- 
librium for diffusion: if condition holds, then the 
time evolution will produce the ensemble average with 
(p(t -> 00)) = and (p 2 (t -> 00)) = (p 2 ) eq - This re- 
sult also suggests that the EH holds, and thus the KT 
holds. Now we may ask in what situation Eq. ([2]) is not 
valid. First, one should note that the long time behavior 
is associated with the small values of z in the Laplace 
transform. Indeed, from the final value theorem [2l| we 
have 



lim R(t) 



lim zR(z) 

Z^O 



(10) 



Therefore, it is only necessary to know R(z). Morgado 
et al. [3] obtained a general relationship between the 
Laplace transform of the memory function Tl(z) and the 
diffusion exponent a: 



lim II(z) w cz c 

z^O 



(11) 



where c is a positive non-dimensional constant. Now, 
using Eq. (fTTj) . in the limit Eq. (fTUj) we have 



lim R(t) = lim (1 + cz c 

t^OO Z-*0 ^ 



(12) 



which is null for all diffusive processes in the range 
< a < 2. In fact, this occurs in equilibrium or near- 
equilibrium states in which the validity of Linear Re- 
sponse Theory holds. On the other hand, this condi- 
tion fails for ballistic motion, a = 2, in which R(t — * 
00) = 1/(1 + c) and the autocorrelation function C p (t) 
will be non-null for long times. In other words, if the 
ballistic system is not initially equilibrated, then it will 
never reach equilibrium and the final result of any mea- 
surement will depend on the initial conditions. In this 
situation, the EH will not be valid; however, once again 
the KT holds since the violation of the EH was due to 
the violation of the irreversibility condition, Eq. ([2]), as 
predicted by Khinchin. The main consequence of the 
violation of this condition is the presence of a residual 
current, Eq. ([8]). However, the effective current can be 
very small compared to (p(0)} and its value, as any other 
measurable property for ballistic diffusion, will depend 
on the value of c. In other words, the system decays to 
a metastable state and remains in it indefinitely, even in 
the absence of an external field. 

5. Time average and the EH. — The time averages of 
correlation functions are crucial for elucidating the prop- 
erties of dynamical processes and play an extremely im- 
portant role in the ergodic theory and, consequently, in 
physics. For diffusive systems governed by the GLE, we 
will show that the condition R(t — » 00) = is sufficient 
for the time average to be equivalent to the ensemble av- 
erage, i.e., for the system to be ergodic. For macroscopic 



3 



systems with a large number of degrees of freedom, the 
effect of past values of the forces usually vanishes for a 
sufficiently large t, and the aforementioned condition is 
quite reasonable. 

Let us consider the time average integral 



L n = lim — 

T^oo T 



Jo 



X{t,t')dt'dt, 



(13) 



For stationary systems, x(*>^') = x(t ~ we arrive 
at 0,3 



lim 

t — >oo 



X (t')dt' + R(t) 



R(t')dt' 



(14) 



Given that R(t) is a real-valued function that converges 
asymptotically to a finite value, since we arc working with 
the velocity autocorrelation, we can use a generalization 
of the final- value theorem for Laplace transforms [211 ], 



lim — 



lim zR(z) 

2^0 



With this, we obtain from Eq. (114 



R(t)dt. 



X(0) + R(t -> oo) 



lim zi?(z) 



(15) 



where Xs is the time independent value, often called static 
susceptibility. On the other hand, taking the Laplace 
transform of Eq. ([7]), we obtain x(z) + zR(z) — Xs- Tak- 
ing the limit z — ► 0, the previous relation becomes 



x(0) + lim zR(z) 



Xs- 



(16) 



Comparing Eq. (fill)) with Eq. (|15p . one should conclude 
that the EH can only be valid if R(t — > oo) = 0, i.e. if 
the irreversibility condition © holds. From Eq. (TTUJ) we 
end up with 



X(0) = Xs 



(17) 



Again this is a consequence of the irreversibility condi- 
tion. Therefore, irreversibility is a necessary and suffi- 
cient condition for the EH to hold in diffusive processes 
described by a GLE. 

6. Simulation. — In order to illustrate the analytical 
results, we have numerically integrated the GLE, Eq. |(SJ), 
to obtain approximations to the probability distribution 
of particle velocities using histograms. We construct the 
memory using 



/ 



n(<) = / p(lu) cos(Lut)duj, 



where 



clujP , for uo < uo s 
g(u>), otherwise. 



(18) 



(19) 





-1.0 -0.5 



0.0 0.5 
P 



FIG. 1: (Color online) Numerical results for the probabil- 
ity distribution function for subdiffusion (top, a — 0.5), nor- 
mal diffusion (middle, a = 1) and superdiffusion (bottom, 
a = 1.5). The time averages (circles) are obtained by follow- 
ing one particle trajectory and calculating the histogram for 
times from t — 100 to t = 5000. For the ensemble averages 
(squares), we calculate the histogram using 5 TO 4 particles, at 
time t — 1000. The continuous line is the Maxwell-Boltzmann 
distribution. Insets: Curves a correspond to the functions t a 
and curves b to the simulated mean square displacements. 



Here, the function g(ui) is arbitrary as long as it is suffi- 
ciently well-behaved and that its integral in the memory 
function converges. If one is interested only in the long 
time behavior t 1/Vs, it can be taken to be 0. With 
this noise density of states, it is possible to simulate many 
diffusive regimes [l5[ . Noise of this form can be obtained 
either by formal methods or empirical data. Using this 
expression in Eq. (|18[) . taking the Laplace transform in 
the limit z — > 0, we have 



U(z) 



z p , for (3 < 1, 
-az m(z), for /3 = 1, 
z, for P > 1. 



(20) 



Consequently, for this type of noise, there is a maximum 
value of a, i.e., a < 2 for any value of (3. It should be 
noted that the case (3 = 1 does not lead to a memory 
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whose Laplace transform is in the form of Eq. (jlip . For 
— 1 < j9 < 1, we obtain a = 1 + 0. For > 1, one has 
a = 2, which shows that ballistic diffusion is a limiting 
case for the GLE with this type of memory. 

In Fig. (TT]) we show the probability distribution func- 
tions obtained for subdiffusion (/? = —0.5), normal diffu- 
sion (/3 = 0), and superdiffusion (f3 — 0.5) for the values 
a = 0.25 and g{w) — 0. We have used uj s — 0.5 for 
all cases except for subdiffusion, which demands a broad 
noise u> s = 2 to reach the stationary state. In all cases, we 
expect that R(t — > oo) = 0, and that the EH will be valid 
even for the subdiffusive (superdiffusive) case, despite the 
fact that W = iW = oo). These relations can be seen 
by considering the limit W = lim z ^o R(z). If the EH 
is valid, the velocity probability distribution will be the 
same for an average over an ensemble of particles and for 
a time average over the trajectory of a single particle for 
long times after the system has reached an equilibrium 
state. Note that despite the presence of large fluctua- 
tions in the time average case due to numerical errors, 
there is a good agreement between the resulting ensem- 
ble and time distributions. The three probability distri- 
butions converge toward the Maxwell-Boltzmann distri- 
bution, which is in accordance with previous analytical 
results [l7| • 

7. Concluding remarks. — In this work, we have shown 
that the KT (proved by Khinchin for normal diffusion) 
holds for all kinds of diffusive processes, which are ergodic 
in the range of exponents < a < 2^ This result may 



have deep consequences in many areas @, IE H, US 11 , 12j . 
Moreover, it could be verified in and applied to experi- 
mental systems, such as the subdiffusive dynamics of the 
distance between an electron transfer donor and accep- 
tor pair within a single protein molecule [22}, which has 
recently been modelled by a GLE (23[. Such a model 
successfully explains the equilibrium fluctuations and its 
broad range of time scales, being in excellent agreement 
with experiments. The KT gives the EH a practical 
character, since it is expressed in terms of response func- 
tions: our results apply for real-valued relaxation func- 
tions R(t); on the other hand, if the relaxation function 
assumes complex values, e.g. conductivity, the final value 
theorem may not be applied. For those systems, the KT 
fails, as proposed in Ref. [l[. 

In principle, it is generally possible to derive a GLE 
for Markovian systems by eliminating variables, whose 
effects are incorporated in the memory kernel and in the 
colored noise [2J]. Altogether, some results obtained for 
the GLE formalism should be valid for diffusion described 
by fractional Fokker-Planck equations, since both for- 
malisms yield similar results |23[ . For nonlinear Hamilto- 
nian maps (l3| there is no general framework to address 
the problem and, in particular, the absence of a coupling 
to a thermal bath (explicit in the GLE) and consequently 
the lack of a detailed balance relation or FDT may re- 
quire a specific analysis of each case. However, since it 



is possible to give a kinetic description of the Hamilto- 
nian dynamics by means of a fractional Fokker-Planck- 
Kolmogorov equation 25J, it is expected that the treat- 
ment of anomalous diffusion in such systems should also 
be possible by the GLE formalism. Further research in 
this direction is needed and will open new perspectives. 
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